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On the corrections to Strong-Stretching Theory for end-confined, charged 
polymers in a uniform electric field 

Gabriele Migliorini 1 

Department of Mathematics, University of Reading, 
Whiteknights, Reading RG6 6 AX, United Kingdom 

We investigate the properties of a system of semi-diluted polymers in the presence of charged groups and 
counter-ions, by means of self-consistent field theory. We study a system of polyelectrolyte chains grafted to 
a similarly, as well as an oppositely charged surface, solving a set of saddle-point equations that couple the 
modified diffusion equation for the polymer partition function to the Poisson-Boltzmann equation describing 
the charge distribution in the system. A numerical study of this set of equations is presented and comparison 
is made with previous studies. We then consider the case of semi-diluted, grafted polymer chains in the 
presence of charge-end-groups. We study the problem with self-consistent field as well as strong-stretching 
theory. We derive the corrections to the Milner-Witten-Cates (MWC) theory for weakly charged chains and 
show that the monomer- density deviates from the parabolic profile expected in the uncharged case. The 
corresponding corrections are shown to be dictated by an Abel-Volterra integral equation of the second kind. 
The validity of our theoretical findings is confirmed comparing the predictions with the results obtained 
within numerical self-consistent field theory. 
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I. INTRODUCTION 
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The stability of a dispersion against coagulation has 
been one of the primary subjects of research in colloid 
science. Colloidal stability can be obtained in polar sol- 
vents by means of ionic molecules. In this case an elec- 
trical double layer forms at the interface between the 
colloidal particles and the solvent, and the dependence 
of its structure on the electrolyte properties has been 
the key ingredient to explain electrostatic-induced col- 
loidal stability in the past 1 . In practice however elec- 
trostatic stability is difficult to obtain, due to extreme 
sensitivity of the process to the electrolyte conditions. A 
viable alternative to electrostatic stabilization is to ab- 
sorb or end-graft non-ionic polymers on the surface of 
colloidal particles 2 -! 3 -. The stabilization mechanism ob- 
tained in this case is steric in nature. The theory of 
stabilization of colloids by adsorbed neutral polymers 4,5 
is one of the first problems where self-consistent field the- 
ory has been applied^ 7 -. Self-consistent field theory, in- 
troduced in the context of neutral polymer systems by 
Edwards^ and Helfand 9 , together with its discrete, lat- 
tice formulation counterpart, has been applied and ex- 
tended to a variety of inhomogeneous systems, includ- 
ing block-copolymer systems^, polymers at interfaces, 
as well as non-ionic polymer brushes^ 2 - - — . The study of 
polymer chains terminally attached to an interface has 
been branching from the colloidal science area to many 
other domains of research. Polymer brushes are a topic 
of interest in surface science, in view of their ability 
to modify surface properties^ 5 - - — . In particular chain- 
end- functionalized polymer brushes, with uncharged 18 , 
as well as charged^— end-groups, have been used to 
design and control surface properties in a variety of dif- 



ferent ways. Theoretical studies of charged-end-group 
brushes^ 2 -, and their response to an electrical external 
field^ 3 -, have been discussed recently. 

Polyelectrolytes are macromolecules containing ioniz- 
able groups. Their tendency to dissociate in solution into 
charged groups and low-molecular-weight counter-ions 
makes polyelectrolytes an interesting problem, as both 
steric interactions of the high-molecular- weight structure 
and electrostatic interactions between counter-ions and 
charged groups, coexist at different length scales. Due to 
the intricate interplay of these different type of interac- 
tions, the properties of polyelectrolyte solutions, as well 
as the behavior of polyelectrolytes near an interface, have 
been considerably more difficult to explain than those 
of non-ionic polymers. Extensions of the lattice self- 
consistent field theory for neutral polymers 10 have been 
discussed 2 -^. The self-consistent theory of inhomogeneous 
polyelectrolyte systems has been formulate d 25 ' 26 , and re- 
cently reconsidered in the context of symmetric di-block 
polyelectrolytes and polyelectrolyte blend s 27 ' 28 . As dis- 
cussed in the paslj 2 ^— , for a system of end-grafted poly- 
electrolytes, the charge distribution is localized at the 
interface and symmetry considerations, as well as con- 
siderations about the double layer that forms at the sur- 
face, simplify their study. We consider both the case of 
end-grafted polyelectrolyte chains and end-grafted poly- 
mer chains with charged-end-groups at a flat interface, by 
means of self-consistent field theory. Related studies have 
been presented in the pas t 22 ' 32- — and a detailed compar- 
ison of the results obtained within self-consistent field 
theory with the predictions of analytic strong-stretching 
theory 3 - 5 - are discussed in this work. The conformational 
properties of the system in the presence of a uniform 
external field^ 3 - are then addressed. The response of a 
polyelectrolyte brush to an external electric field repre- 



sents an important problem with several practical appli- 
cations. Its role in the actuation of nano-size cantilevers, 
to cite one of the many possible applications, has been 
recently reported^. We study a system of end-conhncd 
charged polymer chains and its response to a uniform 
external field. We discuss the conformational changes of 
the monomer-density and counter-ion distribution pro- 
files, for both the case of a similarly and oppositely 
charged grafting surface. We propose a new mathemati- 
cal framework, that generalizes strong-stretching theory 
for charged grafted polymers in the presence of a uniform 
field and we discuss the results of our findings. 



II. THEORY 

We introduce in this section the theoretical methods 
to describe the properties of a polyelectrolyte system of 
chains attached to a flat interface. The system is com- 
posed of n p polymers chains of length N p , grafted to a 
flat, charged substrate of area A with grafting density 
a = n p j A, in the presence of a second surface, at a dis- 
tance d and in the presence of n c counter-ions of valence 
z c in solution. We consider the case of strongly dissoci- 
ating polyelectrolytes, so that the number and position 
of the charges along the polymer chains is fixed and de- 
termined by the ionization degree parameter /, the to- 
tal charge per chain being fN p . The charge distribu- 
tion along the polymer will be referred as z p (s), s being 
the parameter that measures distance along the polymer 
backbone. Both the uniformly charged as well as the 
neutral brush with charged-end-groups are described by 
the present formalism. The ionization degree, polymer 
length and grafting density determine the counter-ion 
concentration, according to the electro-neutrality con- 
dition, n c z c + fN p n p z p = 0. The presence of mono- 
valent salt may easily be included in the notation in- 
troduced below, as well as the presence of other 'small' 
and polymeric charged species, according to the general 
formulation 2 ^ discussed in Appendix [SJ Each polymer 
occupies a volume N p / po and has a natural end-to-end 

1/2 

length of aN p , a being the Kuhn segment. We as- 
sume all charged species, namely counter- and co-ions, 
to be point-like particles, neglecting steric contributions 
due to non-polymeric species. A general formulation of 
the theory, where several species are present and the size 
of the counter-ions and co-ions is explicitly taken into 
account, is discussed in Appendix The system is im- 
mersed in a solvent. As in the case of non-ionic sys- 
tems, the corresponding degrees of freedom can be in- 
tegrated out, so that the effective interaction between 
the monomers is described by the effective, short-ranged 
excluded- volume interaction parameter v. The Bjerrum 
length lg = e 2 /e(r)fcsT, together with the reference den- 
sity po and the excluded-volume parameter v, character- 
izes the properties of the system. 

We consider in this section a single polyelectrolyte 
species in the presence of counter-ions, so that N p = N 



and n p = n. The dimensionless counter-ion concentra- 
tion is defined as 



0c (r) 



Pa 






5(v 



and the monomer concentration is defined as 

n „i 



at r 

k(r) = -£ / S{r-r;(s))ds, 
Po a =i J o 



(1) 



(2) 



where the configuration of the a-th chain, belonging to 
the polymer species p is defined by the vector r" and 
where the position of the a-th counter- ion is given by r" . 
The partition function associated with a system of 
polyelectrolyte chains in the presence of counter-ions, 
tethered to a grafting surface with density a, can be 
converted in a statistical field theory^, according to 
the general formulation of model F and model K, as 
defined in the classical inhomogeneous polymer system 
classification^. After integrating out the solvent degrees 
of freedom, and considering the case of strong polyelec- 
trolytes, so that the polymer charge distribution is inde- 
pendent on the location of each monomer and smeared 
along the chain with a uniform charge ratio /, the effec- 
tive hamiltonian of the system reads 



-^(r) 2 + -i-|V^r)| 2 

V 47UB 



PU[w p ,ip] = | / dv 



—a I dr± In Q p (t±, iw p + iz p ip) — n c In Q c (inp), (3) 



where w p (r) and i/j(t) are the two fields related 
to monomer concentration and electrostatic potential, 
where Q p and Q c are the single-chain and counter-ion 
partition functions, defined and in Appendix [X] and 
where a uniform distribution for the grafting points is 
assumed. The excluded-volume parameter, also intro- 
duced in Appendix \X^ represents the effective interac- 
tion between monomers. As both the electrostatic and 
monomer density fields are uniform in the x and y direc- 
tions, parallel to the grafting surface^, the single-chain 
polymer energy contribution, along the direction perpen- 
dicular to the grafting surface reads 



E[z a -,o,s] 



2a 2 N 



\z a (t)\ 2 +w p (z a (t)) dt. (4) 



The grafting density a can be absorbed introducing 
rescaled polymer and counter-ion concentrations, 



<j)p(z) 



ap 



(j> p {z) 



and 



M z ) = r Ji/o iM 2 ) 



faN 



1/2 



(5) 



(6) 



that are now normalized according to 

I'd 

(f> p (z)dz = aN 1/2 

(f> c (z)dz = aN 1/2 . 

The corresponding field equations are 

w p (z) = Acj) p (z) + N c ip(z) 
w c (z) = -Nif>(z), 



(7) 



(8) 



(9) 



where the reduced interaction parameter is defined as 40 , 

vaN 3 / 2 



A 



ap 



(10) 



and is related to the brush height L according 41 - to 
L/aN 1 / 2 = (4A/-7T 2 ) 1 / 3 , where we introduce the num- 
ber of ions per polymer parameter N c = fN, and where 
all length scales are expressed in units of aN 1 / 2 . 



III. SELF-CONSISTENT FIELD THEORY 

The crux of self-consistent field theory is the one- 
dimensional polymer partition function, 



q(z,z ,s) 



T>z a exp 



E[z a ;0,s] 



x 8(z a { S )-z)8(z a (0)-z ), 



(11) 



for a section of s N segments with its ends at z a (0) = zq 
and z a (s) — z. The polymer partition function satisfies 
the modified diffusion equation 11 



—q{z,zo,s) 



a 2 N d 2 

—^--^2-Wp(z)q(z,z ,s) (12) 
6 az A 



where q(z, Zq,0) — aN 1 ' 2 5(z — zq) and q(0,zo,s) = 0, 
and where the field equations © include a contribution 
from the dimensionless electrostatic potential ip(z), that 
satisfies a non-linear Poisson-Boltzmann equation 



1 d 2 



L B dz 2 
with boundary conditions 

oz 

OZ 



i>(z) = 4>c{z) - (f> p (z), 



(13) 



= -1 



= +1, 



(14) 



that describe an external surface charge, where 

L B = [{ze) 2 /eaN l/2 k b T]/aN c (15) 



is the rescaled Bjerrum length parameter and where the 
Gouy-Chapman length is defined as 



A GC = zfaN/ZL B . 

The counter-ion distribution is given by 

aN 1 / 2 
4>c{z) = — — exp(-w c (z)/N), 



Q, 



where 



Qc 



exp(—w c (z) J N)dz . 



(16) 



(17) 



(18) 



Once the solution q(z, Zo,s) to the modified diffusion 
equation (j!2|) has been found, the concentration profile 
for a chain with its free-end at z = zq can be computed 
according to 



P (z;z Q ) 



1 g(e, z,l -s)q(z,z ,s) 

q(t,za,l) 



ds , 



(19) 



where e is a small finite distance from the grafting surface 
introduced according 39 to the standard formulation of 
model K. The free-end distribution is defined as 



aN 1 / 2 
9( z o) = —p. — exp I 



k B T ' 



where 

fo(zo) 
and where finally 

Q P = [ 

Jo 



exp 



lnq(0,z , 1 

/ fo(zo) 
[ k B T 



dzn 



(20) 

(21) 

(22) 



and the monomer density distribution satisfies 
1 



<fi P (z) 



aN 1 / 2 



g(z )(j)p(z;z )dz . (23) 



The set of self-consistent field equations (jl2 j) - (p!3f) have 
been solved numerically, for different values of the 
four fundamental parameters, namely the natural brush 
height L/aN 1 ' 2 , the number of ions per polymer N c = 
fN, the rescaled Bjerrum Length L B and the relative 
surface charge parameter 



f e = E/zfaN, 



(24) 



where E is the total charge on the grafting surface in 
elementary charge units. The free-energy of the system 
can be computed as 

cd 



F 



1 



k B Tn c aN l / 2 N c JQ 

l-d 



9(zo)[ z ^ +\ng(z )]dz 



k B T 



+ - 



1 



o 



dz(f> c (z) \a.<f> c (z)dz 



1 



A 



2aN l / 2 N r 



2aN i/2 I ^(z)[M z ) ~ M z )]dz, 



4> 2 {z)dz 



(25) 



where the first term is the polymer free-energy and is 
obtained as the sum of the translational entropy of the 
chain ends and the average of the stretching energy 



fe(zg) 



= -hxq(e,zoA) - 



1 



aN 1 / 2 



w(z)<fi(z, z<j)dz, 

(26) 
the second term represents the translational entropy of 
the mobile counter-ions and the third and fourth terms 
the interaction energy, due to excluded-volume interac- 
tions between polymer segments and the electrostatic en- 
ergy related to the charge distribution in the system re- 
spectively. Details on the derivation of equation ([25]) can 
be found in Appendix 151 



IV. STRONG-STRETCHING THEORY 

When the polymer chains are stretched it is possi- 
ble to exploit the analogy between polymer field the- 
ory and quantum mechanics and approximate the single- 
chain partition function in equation (|11[) . defined above 
within the path-integral formalism, by considering the 
dominant contribution of the classical path and the fluc- 
tuations around i t 42 ? 43 . The problem in this limit relates 
to the classical mechanics problem of a particle moving 
down an incline- 4 . In the classical limit, the partition 
function q(0, zq, 1) is dominated by the path z a (s) that 
minimizes the energy in equation Q and an estimate for 
the free-energy fo(zo) of a chain extending to z = z , can 
be obtained accordingly: 



q(0, Z , 1) ex exp ( - E[Z a - l]/k B T] 



(27) 



where, from here on, we will use dimensionless units 
Z = z/aN 1 / 2 . Following the general principle of strong- 
stretching theor y 42 i 44 we write, 

\s\Z, Z ) = N c [^v 2 (Z ) + U(Z) - U(Zo)] , (28) 

where the potential U(Z) is to be identified with the poly- 
mer field w p {z) introduced in equation Q. The chain- 
end tension can be written^ 3 - as 



v(Z ) 



(Jlnff(Zp) 
dZ 



and, in dimensionless units, the two constraints 



1 



and 



o S(Z, Z ) 

D g(z 



-dZ = 1, 



S(Z, Z 



dZ Q = (f)p(Z), 



(29) 



(30) 



(31) 



apply. This yields an expression for the dimension- 
less speed at position Z, S(Z, Z ) and the normaliza- 
tion condition in equation (|30[) . corresponding to the 



isochronicity constraint for the effective particle rolling 
down an incline 4(? . In the absence of excluded-volume 
interaction ;) 35 ' 44 , the normalization condition in equation 
(f3T)| can be written as an Abel-Volterra integral equation 
of the first kind—. Its solution leads to a parabolic form 
for the potential U(Z) 42 ' 46 . Neglecting the term related 
to the end-monomer tension, i.e. setting v(Zq) s=s 0, one 
finds the integral equation 



u 



dtf{t){t~U)- 1 ' 2 = l 1 



(32) 



and its solution, in terms of the Riemann-Liouville frac- 
tional derivative of order 1/2, as discussed in Appendix 
ICl reads 



f(U) 



3dZ 
2dU 



1 d 



U 



it dU ju 

The associated potential form is obtained as 

3^ 2 



dt(t-U)- 1/2 . (33) 



U(Z) 



■(H 2 -Z 2 ), 



(34) 



-1/2, 



and the electrostatic potential 4> a {Z) = N c U(Z) fol- 
lows a parabolic form, according to strong-stretching 
theory 44 . We express distances in dimensionless rescaled 
units Z' = Z/N c ' , where H represents the brush height 
in units of Z . The relative charge ratio per chain reads 

Qi = 1 - v> where 



The rescaled counter-ion distribution is given by 



(35) 



N l J 2 r c (Z) = ^^ exp{^(Z)}, (36) 

and the electrostatic potential outside the brush is writ- 
ten as 



^\Z) = -2 In h 1/2 (Z -H + 1/if 
where the brush height 44 is given by 



(37) 



H + ^ V ^ H2 ^" 2 ^ IT* ) = 7, (38) 



the monomer density distribution, 

N^ P (Z) = N^rAZ) + 1/7, 



(39) 



and the counter-ion distribution outside the brush is 
given by the Gouy-Chapman form, 



4> b c (Z)= l -{Z-H+l/H)- 2 . 

7 

The size of the brush can be expressed as 
(Z) _H 2 ^ H , 



H 



27 



-e s 



(40) 



(41) 



Finally the end-monomer distribution can be obtained, 
as detailed in the Appendix [D] solving 4 ^, a similar Abel- 
Volterra equation of the first kind to equation (|30|) . that 
corresponds to the constraint of equation (f3"Tj). One finds, 



iN'J^iZ) 



^ Z{l + ^f )(r{Z)) 



1/2 



-3tt 



+ ^—H 2 exp(r(Z))erf(^(Z)). (42) 



In the following section we will compare the strong- 
stretching theory predictions with the results obtained 
within self-consistent field theory, where the Poisson- 
Boltzmann and modified diffusion equations p^|) - dT^|) are 
solved explicitly by finite difference methods. We con- 
sider both the case of excluded- volume and electrostatic 
interactions being present. It is important to note at this 
point that typical synthetic polymers considered in ex- 
perimental conditions are relatively insensitive to solvent 
conditions and the long range electrostatic interactions 
will dominate in some cases^. 



V. SELF-CONSISTENT FIELD THEORY-RESULTS 

In the previous section we discussed the main assump- 
tions and results of strong-stretching theory and intro- 
duced the fundamental equations for the monomer den- 
sity and electrostatic potential in a system of uniformly 
charged polyelectrolyte chains, tethered to a flat inter- 
face, in the absence of the external field. The numerical 
method, used to solve the mean-field saddle point equa- 
tions discussed within self-consistent field theory, is intro- 
duced in this section and the results throughout discussed 
and compared with the analytical predictions of strong- 
stretching theory 44 . In order to solve the coupled set of 
saddle-point equations ([T2"T) - (|13|) . we used a real-space, fi- 
nite difference second-order method 4 ^. The modified dif- 
fusion equation (jT2"j) was solved using a discrete Crank- 
Nicholson real-space method, as discussed in detail in 
the pasl4£. The full Poisson-Boltzmann equation (TT5|) 
was solved using an iterative Newton-Raphson method, 
where the partition function for the small species is cal- 
culated at each iteration step. The algorithm we im- 
plemented in order to solve the modified diffusion equa- 
tion (Ti"2"j) includes the Anderson mixing method 48 . We 
observe our algorithm to converge efficiently in a broad 
range of parameter values. It is important to mention 
that we observe the range of convergence, namely the 
highest number of ions per chain the algorithm can reach 
for a given value of the Bjerrum length and brush thick- 
ness parameters, to depend in a non trivial way on the 
Anderson mixing method. Avoiding the Anderson mix- 
ing method and decreasing the mixing parameter /i to 
very small values, typically /i = 0.01, improves the range 
of convergence but increases the time of convergence sig- 
nificantly. A typical threshold value we obtained and 
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FIG. 1. The two upper plots represent the monomer den- 
sity profile for the value of the brush thickness L/aN 1 ' 2 — 2 
and for increasing values of the number of ions per chain 
N c = 10 • ■ • 90 and N c = 30, 60, ■ • • 180 , for the values of 
the Bjerrum Length parameter Lb =0.1 and Lb = 0.001 
respectively. The two lower plots represent the correspond- 
ing counter-ion concentration profiles for N c = 20 • ■ ■ 100, 
iV c = 30, 60, • ■ • , 180. 



report for comparison: N" lax ~ 200, for values of L = 2 
and Lb =0.1 of the brush thickness and Bjerrum length 
parameters respectively. 

A first set of numerical results has been obtained for 
values of the brush thickness parameter L = 2 and Bjer- 
rum length values Lb = 0.1 and Lb = 0.003. The re- 
sults for the monomer density and counter-ion distribu- 
tions, for increasing ionic strength N c = 10, • • • , 100 and 
N c = 30, • • • 180 for two different values of the Bjerrum 
length parameter Lb = 0.1 and Lb = 0.001 are shown in 
Fig [TJ Another set of calculations have been performed 
at the values of L ~ 1.88, Lb — 0.68 to confirm the 
agreement with closely related calculations 33 . We note 
that, when exceeding the values of the number of ions 
per chain N c = 30, considered in that paper as a limit- 
ing threshold, a peak at low values of the distance ap- 
pears, already discussed 4 ^. Our findings resolve the in- 
consistency reported^ 3 , for the monomer density profiles 
at small values of the grafting distance z/aN 1 / 2 ~ 0.05, 
where a peak, as can be seen by close inspection of Fig.Q] 
appears. The presence of this peak has already been 
reported 4 ^ when considering high values of the ion-per- 
chain parameter N c > 50. 

All calculations of this and the following sections were 
performed considering value of the rescaled distance D = 
d/aN 1 / 2 between the grafting surface and the electrode 
to be large enough, so to measure a vanishing counter- 
ion density at the boundaries. In order to obtain such 
condition, at the lowest values of the Bjerrum length pa- 
rameter considered, e.g. Lb = 0.0003, we needed to 
consider values of order D = 10 3 . This was possible as 
the two grids we consider to solve the modified diffusion 
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FIG. 2. Solid lines represent the self-consistent field theory 
results for the monomer concentration profiles for values of the 
brush thickness parameter L/aN 1 ' 2 = 2 and Bjerrum length 
Lb ==0.1, for increasing values of the ion number parameter 
N c = 10, --,90. The dashed lines are the results obtained 
with the numerical strong-stretching approach discussed in 
the text. 



FIG. 3. Brush thickness h/aN 1 ' as a function of the num- 
ber of ions per chain N c , for different values of the Bjerrum 
length Lb = 0.003, 0.01, 0.03, 0.1, neglecting excluded-volume 
interactions. The two dotted lines are the scaling predictions 
corresponding to the osmotic and Pincus regimes^ respec- 
tively. 



equation (IT^|) and equation (ITSl) did not have the same 
size, the second grid exceeding the size of the first, de- 
pending on the Bjerrum length values considered. This 
can be seen looking at the range of distances z/aN 1 / 2 we 
plot in Fig. Q] for different values of the Bjerrum length 
parameter Lb- 

Another method introduced in the past to avoid the 
numerical self-consistent field theory approach is to 
linearize^ the Poisson-Boltzmann equation (flUl) and to 
modify the strong-stretching theory accordingly. As a 
third, intermediate option, we considered explicitly the 
Poisson-Boltzmann equation (TP3")) and combined its nu- 
merical solution with the analytical form of the potential, 
as predicted by strong-stretching theory. We hence as- 
sume a parabolic form for the effective potential w(Z) 
and solve numerically the Poisson-Boltzmann equation 
(|13[) . and an estimate for the monomer concentration 
profile, together with the electrostatic potential, follows 
directly. The procedure is iterated until convergence is 
achieved. A comparison with the results obtained by 
the full solution of equations (|12p and (|13p is shown in 
Fig. [3] Closely related methods and results have also 
been discussed 31 . 

We now consider the case where excluded-volume 
interactions are neglected^. The validity of this 
assumption^ is tested and comparison between polyelec- 
trolyte scaling and strong-stretching theory and the nu- 
merical results obtained within self-consistent field the- 
ory, when solving the saddle point equations (fP2|) - ([T3|) is 
presented. We assume here the relative strength of the 
excluded-volume interactions to be small with respect to 
the average electrostatic interaction strength, so that the 
solvent can be considered as a 6 solvent for uncharged 



chains. This assumption, already discussed and justified 
in detail 35 , will be considered in the results of Figs. [3]- El 
so to compare with scaling and strong-stretching theory. 
In Fig. [3] we present the results for the brush thick- 
ness h/aN 1 ' 2 , as obtained from the first moment of the 
monomer concentration profile, for increasing values of 
the number of ions per chain parameter N c and for van- 
ishing values of L/aN 1 / 2 . The classical scaling predic- 
tions, discussed in the literature 2 ^, are clearly confirmed, 
given in particular the large values of the number of 
ions per chain parameter our algorithm was capable to 
achieve. The data fit the asymptotic behavior expected 
at large values of the number of ions per chain parameter 
N c , as shown by the upper dotted line of Fig. [31 where 
the condition of static equilibration of the osmotic pres- 
sure dictates 2 ^, in the strongly charged regime, the scal- 
ing law h/aN 1 / 2 — N c ' . According to the results pre- 
sented in Fig. UJ an intermediate cross-over regime, the 
so-called Pincus-regime22, is present. The curves, corre- 
sponding to different values of the Bjerrum length param- 
eter Lb — 0.003, 0.01, 0.03, 0.1, cross-over to the interme- 
diate scaling regime, according 29 to the scaling expression 
h/aN 1 ' 2 = N 2 Lb, corresponding to the slope of the sec- 
ond dashed line in Fig. [3J The region for this crossover 
to occur is quite narrow: all curves in Fig. [3] converge to 
the neutral brush thickness 1 — value ho/aN 1 ! 2 ~ 0.54281, 
thus explaining why evidence of the Pincus regime in 
numerical simulations turns out to be an elusive task^S. 
In Fig. [5] we compare the numerical results obtained by 
the numerical self-consistent field theory with the predic- 
tions of strong-stretching theory, as described by equa- 
tion (|39p , for the four values of the strong-stretching pa- 
rameter LbN c ' = 3, 10, 30, 100. At large values of 7 the 
agreement between the numerical results and the theory 
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FIG. 4. Monomer concentration profile as a function of 
rescaled distance from the grafting surface ZjaNj , for dif- 
ferent values of the strong-stretching parameter LbN c — 
3, 10, 30, 100. Solid lines are the results obtained with self- 
consistent field theory. The dashed lines are the predic- 
tions of strong-stretching theory. Agreement between strong- 
stretching and the self-consistent field theory becomes very 
good for large values of the parameter 7. 
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FIG. 5. Brush thickness as obtained from the first mo- 
ment of the monomer concentration profile predicted by self- 
consistent field theory (solid-lines), for four different values 
of the Bjerrum length L B = 0.0003, 0.001, 0.003, 0.1 and the 
prediction of strong-stretching theory in equation (|41[) as a 
function of the strong-stretching parameter 7. The agreement 
becomes very good at large values of 7. A small discrepancy 
though is observed. The reason for this discrepancy is not 
clear. 



becomes very good, except at small values of the rescaled 
distance Z/N c < 0.05, and at the edge of the brush, as 
expected. 

A few conclusive remarks, while closing this section, 
are here in order. As shown in Fig. [51 the fundamental as- 
sumption of strong-stretching theory becomes more accu- 



FIG. 6. The left upper plot represents the end-segment dis- 
tribution g(zo)- The solid curve corresponds to the results 
obtained within self-consistent field theory, at values of the 
strong-stretching parameter LbNJ — 3,10,30,100. The 

dashed line corresponds to the prediction of equation Q42p for 

3/2 
LbN c — 100. The agreement is good for all values of the 

distance from the grafting surface z/aN 1 ' 2 , except at large 
values where the prediction of equation (|42[) diverges. In the 
lower left plot we show the chain-end tension for the four val- 
ues of the strong-stretching parameter above. The upper right 
plot shows the segment profiles (f>(z, Zq) for LbN c = 100 for 
the four values of the end position zo/aN 1 ' 2 — 1, 2, 3, 4. The 
corresponding prediction of strong-stretching theory is repre- 
sented by the dashed lines. The corresponding average poly- 
mer trajectories are shown in the lower right plot together 
with the prediction of strong-stretching theory, shown by the 
dashed line. 



rate for increasing values of the strong-stretching param- 
eter LbN c ' . At values of LbNJ = 100, the chain-end 
tension, as measured within self-consistent field theory, 
approaches zero, for all values of the distance z/aN 1 / 2 
inside the brush. Fig. [5] also shows the segment pro- 
files <P(z,zq), for chains with different end-positions zq 

3/2 
at values of LbN c ' = 100 and the corresponding av- 
erage polymer trajectories z a (s)/z as obtained within 
self-consistent field theory, together with the prediction 
of strong-stretching theory. We finally present the results 
obtained within self-consistent field theory, when consid- 
ering a uniform charge on the grafting surface. We ob- 
tained results for both the case of similarly and oppositely 
charged surface. The monomer-density and counter-ion 
distribution profiles are shown in Fig. [7] and Fig. |5J 

The results were obtained for both positive and neg- 
ative values of the relative surface charge parameter 
/ e , for a brush thickness parameter value L/aN 1 / 2 — 
2. The boundary condition in equation ([TO) has been 
changed accordingly. As expected, for the case of a 
similarly charged surface, the polymer brush swells, and 
the monomer-density profile obtained at small values of 
the relative surface charge parameter changes its shape 
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FIG. 7. Monomer concentration <j> v (z) and counter-ion distri- 
bution profiles 4>c(z), for brush thickness values of L/aN 1 ' 2 — 
2 and for two values of the Bjerrum length parameter Lb — 
0.1 and Lb — 0.001, as a function of the distance from 
the grafting surface, for positive vales of the relative sur- 
face charge parameter / e = 0.5, for increasing values of the 
number of ions per chain parameter N c — 10, • • • , 40 and 
N c — 30, 60, ■ • ■ , 150 respectively. 
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FIG. 8. Monomer concentration <j> P (z) and counter-ion distri- 
bution profiles 4> c {z), for brush thickness values of L/aN 1 ' 2 — 
2 and for two values of the Bjerrum length parameter Lb — 
0.1 and Lb — 0.001, as a function of the distance from 
the grafting surface, for negative vales of the relative sur- 
face charge parameter / e = —0.5 and for increasing values of 
the number of ions per chain parameter N c = 10, • • • , 60 and 
iV c = 30, 60, • • • , 240 respectively. 



and deviates from the parabolic form. Similarly, when 
considering an oppositely charged grafting surface, as in 
Fig. HI the brush contracts and the counter-ion distri- 
bution changes self-consistently. Related results, for the 
case of a similarly charged grafting surface, have been 
reported in the past 3 ^. 



VI. CHARGED-END-GROUP POLYMER BRUSHES 

In the previous section we presented numerical self- 
consistent field theory results for uniformly charged, end- 
grafted polyelectrolyte systems, both in the presence and 
in the absence of excluded-volume interactions and dis- 
cussed the effects of a uniform electric field. We discussed 
the results for the chain-end, counter-ion and monomer- 
density distributions, for different values of the relative 
surface charge parameter / e . The analysis suggested that 
the parabolic form of the monomer density profile is dras- 
tically affected by the electric field, both in the similarly 
and the oppositely charged case. 

We consider in this section the case of a polymer brush 
in solution, characterized by charged end-group s 22 ' 23 in 
the presence of counter-ions and excluded-volume inter- 
actions between the monomers. Rather than considering 
a uniformly charged polyelectrolyte, as discussed in sec- 
tion|Vl we study in this section the case of grafted neutral 
chains, with a functional charged end-group, represented 
in our model by the terminal, non-grafted \N monomers 
along the chain. As discussed above, we consider the case 
of strong polyelectrolytic functional groups, so that the 
charges are assumed to stay bounded to the chain-ends. 
As for the case of uniformly charged polyelectrolytes, the 
N c ions per chain are uniformly smeared along the ter- 
minal, non-grafted section of the chain, with a rescaled 
charge density f/x so that, consistently with the nota- 
tion of Appendix fAl we consider the charge distribution 
along the chains to be given by z p (s) — //x@( s ~x)- We 
solve the saddle-point equations (fl2")) - ([lU)) that now in- 
volve both the charged and uncharged monomer densities 
and the effective field w p (z) is also defined accordingly. 
In Fig. [H] we show the total (charged and uncharged) 
monomer density distribution of a brush of thickness 
L/aN 1 / 2 — 4, for increasing values of the number of ions 
per chain-end parameter N c — 1, 10, 20. For small values 
of the ions per chain parameter, namely N c — 1, the elec- 
trostatic interactions between the chain-ends are negligi- 
ble, as can be seen in Fig.[S]and the monomer distribution 
is close to the parabolic prediction of strong-stretching 
theory, the brush thickness parameter L/aN 1 / 2 being rel- 
atively large, as expected. Deviations from the parabolic 
profiles can be observed and are due to the inter-chain- 
end electrostatic interactions at higher values of the ion- 
per-chain-end parameter N c . The calculation has been 
performed for the specific value of x — 1/30, so that 
the size of the polymer chain is about thirty times the 
size of the chain-end group. We performed several test 
runs and checked that for smaller sizes of the chain-end 
group the monomer density profiles did not change sig- 
nificantly. For increasing values of N c the brush swells, 
due to the repulsive electrostatic chain-end interactions 
and the counter-ion distribution changes accordingly. We 
now discuss the response of the system to an electric field. 
In our analysis we consider small values of the number 
of ions per chain-end parameter, namely N c — 1, so to 
minimize the interaction between the chain-ends. 
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FIG. 9. Monomer density profiles for a brush with thickness 
parameter L/aN 1 ' 2 = 4, Bjerrum length parameter value 
Lb — 0.1 for increasing values of the number of ion-per-chain- 
end parameter N c = 1,10,20. The dashed line represents 
the prediction of strong-stretching theory, in the absence of 
charged-end-groups. 



FIG. 10. Chain-end distribution g(zo) for a brush with 
thickness parameter L/aN 1 ' 2 = 4, Bjerrum length param- 
eter value Lb =0.1 for increasing values of the number of 
ion-per-chain-end parameter N c — 1, 10,20. The dashed line 
represents the prediction of strong-stretching theory, in the 
absence of charged-end-groups. 



As discussed in section HVl we consider the brush to be 
grafted to a fiat surface that acts as an electrode, while 
the second electrode is kept at a distance D ^ L. In 
order to study the response of the system to the uni- 
form electric field, we solved the saddle point equations 
(H2J) and (fT3"|) changing the boundary conditions for the 
electrostatic potential as well as modifying the effective 
field accordingly. Fig. [11] shows the monomer density 
profiles for the brush of thickness L/aN 1 / 2 = 4, Bjer- 
rum length Lb = 0.1 and number of ions per chain-end 
parameter N c = 1 for both positive and negative val- 
ues of the surface charge f e Ls = ±2. The neutral sur- 
face charge case f e — is also shown. As expected, we 
observe deviations from the parabolic profile for finite 
values of the electric field strength. Differently, for van- 
ishing values of the surface charge f e — 0, the profile is 
very close to the strong-stretching parabolic prediction, 
the value of the brush thickness L/aN 1 / 2 = 4 being rel- 
atively large^. In the following section we discuss how 
to generalize strong-stretching theory for the case of a 
uniform electric field, assuming the relative charge at the 
chain-ends to be small. 



VII. STRONG-STRETCHING THEORY OF 
CHARGE-END-GROUP POLYMER BRUSHES 

The mathematical foundations of strong-stretching 
theory stem from the fact that the isochronicity con- 
straint in equation (|5U| may be regarded as an Abel- 
Volterra integral equation of the first kind. This deter- 
mines uniquely the form of the potential^. To deter- 
mine how a uniform tension at the chain-ends, induced 
by the electric field, affects the potential form U(Z), 



we evaluated the corresponding corrections, solving the 
isochronicity constraint condition in equation (J30]). The 
presence of the electric field induces a chain-end uni- 
form tension and converts the associated integral equa- 
tion from an Abel-Volterra equation of the first kind to 
an Abel-Volterra equation of the second kind, as also dis- 
cussed in detail in Appendix [Cj This leads to evaluate 
the main corrections to the parabolic potential, due to 
the presence of the field, as we now derive. Let us con- 
sider directly the constraint in equation (;3()}, that can be 
written, following the few steps in Appendix [C] 



u a 



u„ 



dtil-v)- 1 ' 2 = / dtf{t)[ir-a{t- V )- 1/2 }, (43) 



where f(U) — \ \^m is related to the first derivative of 

the inverse function Z(U). We then find corrections to 
the strong-stretching results. Namely, after a few steps, 
we obtain the following Abel-Volterra integral equation 
of the second kind: 

df(v) *" ff s d 



F{rf) 



g(v) — 

a 



dtg(r))(t - vr 1/2 , (44) 



where g(jj) = — 1/a, and where a — v6Vb measures the 
uniform end-monomer force induced by the electric field 
at the chain ends. The solution to this Abel-Volterra 
integral equation of the second kind can be obtained as 
follows: 



/(»?) 



dZ 

dij 



1 



o 



r(U -n) 



Jit \Z or 



o-r?) , (45) 
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T(a, x) being the incomplete gamma function of order a. 
The above expression can be integrated and the corre- 
sponding equation for the potential reads 



The expression above reduces to the parabolic strong- 
stretching potential when the chain-end tension vanishes, 
as expected. From the potential in equation pp| , we 
determined the monomer-density and chain-end distri- 
butions by means of the field equation ([5]) , in the limit of 
weakly charged end-groups. This assumption has been 
justified and discussed in Fig. [TUl A similar analysis, 
detailed in Appendix [C] leads to an explicit expression 
for the chain-end distribution g(Zo), for small values of 
the electric field strength a. We find, after a few steps 
detailed in Appendix [C] 



dZ 2 ]/ 2 a 



arj 



V 



r 3/2 



Equation (|47p reduces to the well known expression for 
the chain-end distribution^, when no tension is present 
at the chain-ends. The value of the chain-end prefactor 
S ~ 0.263 has been obtained comparing the results of 
numerical self-consistent field theory, as in the inset of 
Fig. [T51 and the results obtained from equation PS]) for 
the brush height h/aN 1 / 2 . It would be interesting to 
see if the prefactor can be computed explicitly, and if it 
depends on the brush thickness parameter L/aN 1 / 2 . 

Equations (|I6"|) and ((17)) generalize the MWC strong- 
stretching theory to problem of charge-end functional- 
ized brushes in the presence of a small, uniform external 
electric field. A comparison with the results of numer- 
ical self-consistent field theory is shown in Fig. [TTJ and 

Fig. m 

We conclude observing that the monomer-density pro- 
file is affected by the uniform electric field in a different 
way that what seen for uniformly charged polyclcctrolytc 
brushes. Note however that, when comparing the re- 
sults of this section with the monomer-density profiles of 
Fig. [7] and Fig. [5] one should keep in mind that the re- 
sults discussed for the polyelectrolyte problem were ob- 
tained at relatively large values of the number of ions 
per chain parameter N c , the electrostatic interactions be- 
tween monomers being relatively strong. It is interesting 
to observe that for a uniformly charged polyelectrolyte 
system, at low values of the number of ions per chain pa- 
rameter, e.g. N c = 1, the effect of the uniform electric 
field is very similar to what observed for charged-end- 
group chains and the monomer-depletion effect close to 
the surface, seen in the analysis of Fig. [Jj for a similarly 
charged surface, is not present, as in the case of grafted 
charge-end-group polymer chains discussed in this sec- 
tion. This suggests that the strong-stretching theory cor- 
rections discussed in this work might apply reasonably 
well for the problem of weakly charged polyelectrolyte 
chains in a uniform electric field. 
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FIG. 11. Monomer density profile 4>p(z) for a brush with 
thickness parameter L/aN 1 ' 2 — 4, Bjerrum length parameter 
value Lb =0.1, for values of the number of ion-per-chain- 
end parameter N c = 1, for neutral surface charge f e — 0, as 
well as similarly and oppositely charged grafting surface,with 
rescaled surface charge parameter values feLs = ±2.0. The 
dotted line represents the prediction of strong-stretching the- 
ory, in the absence of charged-end-groups. The dashed lines 
represent the prediction of the modified strong-stretching the- 
ory (|46|) . for values of the parameter a — 8J £ Lb- The 
inset shows the monomer-density at the reference distance 
z m = 0.378, as measured within self-consistent field theory, 
shown by the solid line, and as predicted by equation (|46[) . for 
different values of the chain-end tension parameter a, shown 
by the star symbols. 



VIM. CONCLUSIONS 

In this work we studied the properties of charged poly- 
mers grafted to a flat interface, in the presence of an 
external electrical field. Different charge configurations 
have been discussed. We firstly discussed the properties 
of polyelectrolyte brushes, where the charges are uni- 
formly distributed along the chains, while counter-ions 
are in solution and free to diffuse between the graft- 
ing surface and the second, oppositely charged electrode. 
The monomer-density and counter-ion distributions have 
been obtained by numerical self-consistent field the- 
ory and comparison has been made with the analytical 
strong-stretching theory predictions. The effect of a uni- 
form electric field, both in the oppositely and similarly 
charged case, has been studied. We considered charge- 
end-functionalized, grafted polymer chain systems. As 
for the case of polyelectrolyte brushes, we studied the 
conformational and electrostatic properties of the brush 
by means of numerical self-consistent field theory. In the 
case of weakly charged, end-group-functionalized poly- 
mers we have been able to generalize the predictions of 
strong-stretching theory to the case of a uniform electric 
field. We have shown how the corrections to the MWC 
theory for the monomer-density and chain-end distribu- 
tions can be derived in closed form. The presence of an 
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FIG. 12. Chain-end distribution g(zo) for a brush with thick- 
ness parameter L/aN 1 ' 2 = 4, Bjerrum length parameter 
value Lb =0.1 for values ol the number ol ion-per-chain- 
end parameter N c — 1, for neutral surface charge f e = 0, as 
well as similarly and oppositely charged grafting surface,with 
rescaled surface charge parameter values / e I/B = ±2.0. The 
dotted line represents the prediction of strong-stretching the- 
ory, in the absence of charged-end-groups. The dashed lines 
represent the prediction of the modified strong-stretching the- 
ory in equation ()47[) , for different values of the chain-end ten- 
sion parameter a = Sf e L,B- The inset shows the size of the 
brush h/aN 1 ' 2 , as predicted from numerical self-consistent 
field theory, shown by the solid line, and as obtained accord- 
ing to the modified strong-stretching theory in equation (|46|) . 
shown by the dotted line. The dashed line shows the results of 
the modified strong-stretching theory, shifted by the zero-field 
value, obtained subtracting the brush size as obtained from 
self-consistent field theory and the MWC theory prediction. 



external field converts the corresponding Abel-Volterra 
equation from an integral equation of the first kind to 
an integral equation of the second kind. Deviations from 
the parabolic profile have been obtained. A comparison 
of the theory with numerical self-consistent field theory 
has been performed and agreement has been observed. 
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Appendix A: Partition function and Statistical field theory 

In this appendix we review the basic formalism for a 
multicomponent polyelectrolyte system in solution. The 
saddle point equations (fT2|) and (fT3|) for a system of a 
single species of polyelectrolyte chains in the presence 
of counter-ions, discussed and studied in this paper are 
derived in detail in this Appendix and some simplifica- 
tions, compared to a closely related derivation^ 5 -, are in- 
troduced. The average canonical partition function of an 
inhomogeneous multicomponent polyelectrolyte system 
is writte n 25 ! 38 as a functional integral over the cartesian 
position of the nt small species, namely solvent molecules 
and counter-ions, over the conformation of charged poly- 
mer chains, belonging to n p polymer distinct species, and 
finally over the charge distribution in the system, given 
the dimensionless polymer and small species concentra- 
tions in equations ([1]) and (|2j): 

n v 

< z >=n[n e p(H(s)})}z({ Z ;(s)}), (ai) 



where 



*=n(S)/n(n^)/n(n*y 

j 3 ' J p i J k i 

exp (-^P 2 oJ2J dvdv '^ ^ V n' ( r " r ')h' ( r ')) 
jj' 

x <*( e *i N i n i) n s {J dr ^ {r) _ njNj ) 

xII^E^w- 1 )- (A2) 



In the expression above, the chain conformation path in- 
tegral measure Vr p — dr p P/v (r p ), 



Pn p (tD = exp f - — j / ds[r l p (s) 



2a 2 



N„ 



(A3) 



is the Wiener measure, where N p indicates the length 
of each polymer species, a is the Kuhn segment, Q — 
exp~' 3p ° w '"( ) /A 3 is the reference chemical potential 
for species j, that includes the contact interaction en- 
ergy and depends^ 5 - on the thermal wave-length Xtj — 
h/ yjlixm jksT. The canonical partition function has 
to be averaged over the charge distribution P({z p (s)}), 
characterizing each polyelectrolyte present in solution in- 
dependently. The interaction energy term includes both 
short-ranged contact and long-ranged electrostatic inter- 
actions, 



Vy(r - r') = Wij(r - r') + £y (r - r') 



WijSir - r') + 



fc- A/i Ai -j 



e(r)|r — r'| 



(A4) 



where Zj is the valence of each species and should be con- 
sidered as a function of s for any j that corresponds to a 
polymer species, where Vj = p /p° measures the different 
size of each molecule with respect to the reference density 
p , Wij is the contact energy between species i and j, e 
is the elementary charge unit, e(f) the dielectric function 
and j3 = 1/ksT is the inverse temperature in Boltzmann 
units. The first two constraints in the above expression 
for the partition function relates to charge and particle 
conservation, while the last constraint refers to the in- 
comprcssibility condition. We introduce the set of fields 
Wj(r), (J)j(t) and express the canonical partition function 
in terms of a statistical field theory 3 ^, using the particle- 
to-field transformation, 



(Z) = J HjPujVh] exp [ - PHiiwj}, {&})" 

3 



where Stirling's formula has been used and where the 
canonical partition function is expressed in terms of the 
single species partition functions Qj , as an integral with 
respect to density and auxiliary fields of an effective 
Hamiltonian or 'action', 



j 
+ ^E/drdrV J (r)V^^(r') 



-Po 



iV, Q 



^ N- A,. FU Z-^ N, 



tjV^QNj 



Ni 0. 



Nj po 



(A6) 



where cf>j — NjTij/p V enforces the particle conserva- 
tion constraint, where V is the volume, where the single 
species partition functions are: 



Qk = v / dr exp ( ~ / drkuJk ^} 

Q p = i|w( S ) J2 p ({ z p( s )}) 

{z p (s)} 

{-ij ds[u p (r p (s))j, 



x exp' 



(A7) 
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and where finally we have considered the same charge 
distribution within each species being 



where, 



dsz p (s)P(z p ) = z p f, z k = z k . 



(A8) 



Xij 



(3po 
2 

]Po_ 

2 



. " N 



2Wt- 



j ^ Po 
i^Wu + ^Wj/, . (A13) 



Both the smeared and the annealed charge distribution 
case are included in the present formulation, while the 
quenched distribution problem requires a different for- 
mulation. The saddle point equations corresponding to 
equation (|A6[) , up to a redefinition of the fields Uj are, 

-ioji-^^J dr'Vij(r - r')&(r')dr' + , 7 (r)^ : = 



The saddle point equations for the solvent degrees 
of freedom can be integrated out. In the semi-dilute 
regime 2 -, the saddle point equation for the partition func- 
tion Q s defined in equation (|A10[) can be expanded to 
obtain, 



i(j)i 



fry sq 3 






(A14) 



(A9) 



where rj (r) is the Lagrange multiplier corresponding to 
the incompressibility constraint and where charge conser- 
vation applies: J2j ZjNjUj = 0. The second set of saddle 
point equations, for the polymeric and small species, ex- 
plicitly reads: 



P = 7T / dsq p (r,s)ql(r,s) 
«p Jo 

4>k = 7^-exp{-iw fc (r)} 



(A10) 



where the polymer propagator q p (r, s) is the solution of 
the modified diffusion equation, 

^ / n a 2 N p 2 . 

— q p (r, r , s) = — — V q p {v, r , s) - iuj p (r 7 s)q p {r, r , s), 



g p (r,r o ,l) = (a 2 A0^<5(r-r o ) 



(All) 



Integration of the Solvent Degrees of Freedom 

It is useful to rewrite the effective hamiltonian in equa- 
tion (|A6[) and the corresponding saddle point equations 
for the case where a single solvent species is present^. 
One finds, 

j 

3J' 

+ P4j2[drdr><t> j (r)S jj ,<i> j ,(r>) 
jf 

-^E^r^-^ + ^E^^' ( A12 ) 

3 iVj ^ 3 



where the primed sum is over all species but the solvent 
species. The Lagrange multiplier rj(v) can be written as 

Vi^Vs = Y^ \ V 3 ~ Xja)M T ) ( A15 ) 

so that a new set of reduced saddle point equations, with 
i y^ s, holds 

E ««& w + ^- E / dr '^ ( r ' - r )^ « 

(A16) 



J#s 



«</>; 






where «y = (-^- — Xij)> where the index j is intended 
to run over all species excluding the solvent and where 
finally the Flory parameter between species i and j is 
defined in equation (JA13I) . 



Appendix B: Derivation of Model Details 

We consider explicitly the case where we can neglect 
the excluded-volume interactions of the small species. 
For two distinct indices for the polymeric and small 
species, considering the dielectric constant to be uniform 
across the system and considering the case of a smeared 
polyelectrolyte charged distribution, with charged ion 
fraction parameter /, we find: 



/3e 2 
w*>k(r) = po—z c 



dv- 



Y.i*i<t>i(*') 



e\v — r 

Zr) 



ioj p (r) = V" v pp ,(j) p ,(v) + i—fuJk 



Zk 



(Bl) 



where ioj c is related to the dimensionless electrostatic po- 
tential ip(r) according to tp(r) — ioj k (r)z k . Let us con- 
sider the case of a single polymer species, and let us con- 
sider a single species of point particles: in this case the 
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above saddle point equations reduce to 

iV 2 w c (r) = -4-kIbPqZc ^ Zjfc (r) 



iwp(r) = u^p(r) + i — fuj c (r). 



(B2) 



Appendix C: On Abel-Volterra equations of the second 
kind 



In this section we give a detailed derivation of the inte- 
gral Abel-Volterra equations of the first and second kind 
discussed in section [VlJ Consider the integral equation 



The corresponding mean-field free-energy can be com- 
puted according to equation (JA6I) . The canonical Model 
F field theory partition function^ reads: 

(Z) = Z DojpDuJc exp { - pU[u p , uj c }} (B3) 



where 



PU[uj p ,uj c ]/pq = — / dv{io p -fco c ) 

zv J z c 



V 



r 



and where the partition function, 

Z = exp{0 p ln0p + <j) c liL<j) c + fi pV(j> p + iiQ C V(i>c}, (B5) 
also written as 






tp.Jl C - 



Tin •Tic 



(B6) 



is the ideal non-interacting partition function of a mix- 
ture of n p charged polymers and n c counter-ions in solu- 
tion. Changing variables we write, 

f 1 1 

pH[u,w c ] = / dr—cj 2 + — — -|Vw n |2 



2v 8irl B z 2 ' 

n p In Q p [iu + i — fuj c ] - n c ln Q c [iuj c ] . (B7) 



The saddle point free-energy is 



F 



k B T 
~2v 



/3H[-iw,-itp]/p 



8irlr 



dr\Vip(r)\ 2 



— / drw 2 (r) - n p ln Q p [w p + z p if\ - n c ln Q c , (B8) 



where v is the excluded-volume parameter and where 
Ib = e 2 /k B Te is the Bjerrum length. Note that the elec- 
trostatic interaction energy term can be rewritten, using 
equation (|B2|) . 



K 1 f 1 C 

- n p ln Qp [w p + Zpfip] - n c ln Q c [ip] . (B9) 



Equation (|B9[) . after a few steps can be shown to reduce 
to the free-energy expression discussed in equation (|25[) . 



Ua 



We write! 3 ., 



dtf(t){t-U)- 1/2 =g{U). (CI) 



dUogiUoWo-ri)- 1 / 2 



fU r U 

/ dU (U -v)- 1/2 dtf(t)(t- U)-V 2 . (C2) 
Using Fubini's theorem we find, 



dUg(U)(U- V )-^ 2 



Uo 



dtf(t) / dU{U-i])- 1/2 {t-U)- 1/2 . (C3) 



According to the assumption of strong-stretching theory, 
one neglects the chain end tension Vo ~ 0; the parabolic 
form in equation (|33p follows considering the constraint 
in equation (|30p . corresponding to g(U) — 1 in equation 
(ICTl) . This yields 



/W = -(t/o-r7)- 1/2 , 



(C4) 



that can be integrated to obtain the parabolic form in 
equation ([34]) . Similarly, an expression for the end- 
monomer distribution can be obtained considering 



H 



dZ g(Z ) (Z 2 - Z 2 )- 1 ' 2 = ~<j> p (Z) (C5) 



and inverting this expressio n 13 ! 44 by fractional differen- 
tiation, 

1 d r Uo 
fW = —4f dtg(t){t-U)- l '\ (C6) 



one finds equation (J42J). As discussed in section[VTJ let us 
consider now the case of a uniform tension at the chain 
ends. The integral form in equation (|C3[) can be written 
as 



Uo 



dU{U - n) 



_^-V2 = 



Uo 



dtf{t){\ 



(t-v)-v 2 
(t-v) + v 2 



)}• (C7) 
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Expanding to leading order yields an Abel Integral equa- 
tion of the second kind, 



and the solution to the Abel-Volterra integral equation 
of the second kind (IC9I) has the form 



dU(U-T])- 1/2 = j 



dtf{t){*-a{t-ri)- 1 ' 2 ), 



where a = v6| Vq| measures the uniform end- monomer 
tension induced by the electric field, as discussed in sec- 
tion Ell We find: 



m 



a d r Ua 
ndr] J n 

l_d_ 
7T dr\ j v 



U 



dtf(t)(t- V )- 1/2 

dtit-r!)- 1 ' 2 . 



This can be written as 



/(»?) 



1 d 

7T dr\ j v 



U 



and 



f.'.i 



/(*)(* -ri)- l l 2 dt = af{rj)-l. 



(C9) 



(l + af(t))(t-v)~ 1/2 dt, (CIO) 



(Cll) 



fin) = F(r,) 



TV 

772 



Uo 



dte~ 



r(Uo-t) 



F(t) (C14) 



(C8) and after a few steps we find 



m 



-(u -n) 



rA 



1 7T 

2' a 



2 (Uo-v)), (C15) 



where T(a,x) is the incomplete gamma function of or- 
der 1/2 and where the potential form can be obtained 
integrating equation (|C15[) . We note that the function 



/([/) = J ^dZ/dU can be written as 

f{ri) = ^ r u(\,\,^ i {U -ri)), (C16) 

where ^(^,^,-^(Uo — f])) i s the confluent hyper- 
geometric Tricomi function and where one recovers the 
standard strong-stretching theory prediction to first or- 
der expansion. It is easily checked that the first order 
in the above mentioned expansion corresponds to the 
strong-stretching theory result: 



m 



-(Uo-vr 1/2 , 



that corresponds to the parabolic form (|34 
(|C15[) can be integrated, and one finds 



I- 



t(Uo-U) x 

r(-,i)e'dt, 



(C17) 
Equation 

(C18) 



The presence of the external field converts the Abel- 
Volterra integral equation of the first kind (IC1[) . as usu- 
ally discussed within strong-stretching theory, to an inte- 
gral equation of the second kind (|C9[) . The integral equa- 
tion (JC9I) can be solved and a new form of the potential 
follows, as we will show in the rest of this Appendix. To 
leading order in the Taylor expansion, we find 



dm 

dij 



7T , , d 

or dg 



(C12) 



F(v) = g(v) - - 

a 



Uo 



dtg(r))(t-r)) 



-1/2 



The isochronicity constraint in equation (1311) requires the 
function g(rf) to be constant, namely g(rj) = —1/a, so 
that 



F( V ) = -- + ^(U -v) 1/2 , 
a or 



(C13) 



and finally 



|ttZ = 2(U - U) 1 ' 2 



a i= ei 

\/7T 



-(Uo-U) 



r(l^(u -u)). 



(C19) 



The above expression has been used to obtain the poten- 
tial, for both negative and positive values of the chain-end 
force, shown in Fig. [TT] and discussed in section |VT1 For 
the chain end distribution we find 



ds(y) 
drj 

G(v) 



a z dg 

1 
a 3a 2 



(C20) 



1 3/2 



where s(g) = g(g)^jj-. The expression for the chain-end 
distribution in equation (|47|) can be obtained solving the 
above second kind Abel-Volterra integral equation as well 
as using the result of equation (|C15[) above. 



